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Abstract 

The paper is devoted to the analysis of Fermi acceleration of pro- 
tons in GRBs and its neutrino signature. We have compared the 
consequences of Bohm scaling and those of a Kolmogorov scaling, the 
latter being more reliable. The predictions about the energy limitation 
of UHE-protons by the various losses and the neutrino emissions turn 
out to be very sensitive to these scalings. We consider Kolmogorov 
scaling as the most realistic and predict a reasonable pp-neutrino emis- 
sion around 100 GeV at the end of the radiative stage of the fireball 
expansion, for a large number of GRBs pending on their baryonic 
load. A second p7-neutrino emission is expected with the accelera- 
tion of protons in the radiation free stage, but with a synchrotron loss 
limitation immediately followed by a severe expansion loss limitation. 
According to the Kolmogorov scaling, the protons could not reach the 
UHE-range. Anyway the large possibility of a two component neutrino 
emission would be an interesting clue of cosmic ray physics. 

Keywords: GRBs, Cosmic Rays, Neutrinos, Fermi processes 



1 Introduction 

Gamma Ray Bursts (hereafter GRBs) are unique high energy phenom- 
ena in astrophysics because of their possibility to manifest all the interest- 
ing " astroparticles" processes, such as generation of high energy gamma 
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rays, ultra high energy cosmi c rays 
tional waves (see 



. . ^ high energy neutrinos and gravita- 

, Dermei] (|2nni[ ) and iMeszarol (Hool)). The "fireball" model 

( R,ees fc Meszaroslll992h has been successfu in explaining the afte rglow stage, 

and it s more elaborated form with the addition of internal shocks ( Rees fc Meszaro'i 
( 1994 ). jPaczvhski &: Xu (1994)) during the "free" exp ansion has been suc- 
cessfu l in explaining the light cur ve and gamma spectra (jPishman &: Meegan 
ISeloborodov et alT()200(J )). Because of the highly relativistic dynam- 
ics necessarly involved in the GRB phenomenon, the generation of very high 
energy particles is expected through strong shocks and strong magnetic per- 
turbations. Indeed GRB population could be a main source of UH E-cosmic 
rays which coul d be generated by the external shocks ( Vietri 1995[ l or by the 
internal shocks ( Waxman 1995| ). and a fiux of neutrinos produced by the col- 
lisions of the UHE-cosmic r ays with the GRB gamma photon is reasonably 
expected fjWaxman &: BahcalL1999l . 

In this paper, our intend is to look at the sensitivity of the predictions 
of cosmic ray and neutrino generation to the description of the Fermi accel- 
eration process together with opacity effects for protons. We have analysed 
the consequences of two assumptions, one is the so called "Bohm scaling" 
assumption, the other is what we called "Kolmogorov scaling". The former 
consists in the statement that the Fermi acceleration time is proportional to 
the Larmor time of the accelerated particle, with a constant proportionality 
factor larger than unity (often chosen between 1 and 10 in the litterature). 
The latter consists in taking into account that the ratio between these two 
characteristic times depends on the rigidity of the particle through a law 
governed by the t urbulence spectrum, as is confirmed by numerical works 
fICasse et al.ll200lh , where no Bohm scaling has been found. In astrophys- 
ical media such as the solar wind, the interstellar medium, the turbulence 
spectrum is likely consistent with the Kolmogorov law. 

The Bohm scaling is very convenient to make a first investigation of the 
high energy physics performances of the objects, and our prejudices are often 
grounded on this. However, in this paper, we will show that this scaling leads 
to unrealistic results. The analysis based on the Kolmogorov scaling is more 
reliable, and it turns out that it leads to very different conclusions, as we 
will present with details in the paper. The analysis is developed both in the 
radiation free stage and in the radiative stage. We especially emphasise the 
regime when the fireball is opaque to pp-collisions during the beginning of 
the internal shocks stage because of its baryonic load, which should occur 
for a large fraction of GRBs population. We make unusual predictions about 
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the performance of cosmic ray generation, and about neutrino emission that 
could have two components, namely a non-thermal pp-neutrino emission and 
a p7-neutrino emission. 

The paper is organised as follows. In section 2, we briefly present the dy- 
namical description we need for the estimation of the opacity effects, not only 
for photons but also for protons. In section 3, we analyse the consequence of 
the fast acceleration regime governed by the Bohm scaling. In section 4, we 
present the properties of the more progressive acceleration process governed 
by the Kolmogorov scaling and predict a reasonable non-thermal pp-neutrino 
emission with its spectrum. We investigate, more briefly in section 5, the ra- 
diation free stage, where UHE-cosmic rays are expected to be accelerated 
and p7-neutrinos generated. We end the paper with a discussion that sum- 
maries our conclusions about cosmic ray and neutrino generations and the 
sensitivity of the predictions to the acceleration model. 



2 Preliminary considerations 

2.1 Dynamics of the fireball 

In this subsection, we summarize all th e results we need for this paper 
that describe the expansion of the fireball fiMeszaros et al. 1993| ). The wind 



flow is considered to be a set of discrete shells which are successively emitted 
with an energy Eg = E/Ng, where Ng is the total number of shells. The 
duration, ty^, of this wind flow provides with an interval of shell number, 
namely, 1 < A^s < ct^/ro where tq is the size of the central object. 

At the very beginning of the expansion of a shell, the pressure is supposed 
to be dominated by the radiative pressure. The temperature, T, of the plasma 
which is mainly composed by electron-positron pairs is equal to the photon 
temperature. Considering a shell is initially spherical with a radius, ro, and 
with an energy emitted in 7-rays equal to Eg, we have 

(4/3)vrr3-" ' 
where a is the Stefan constant. Thus, the plasma temperature is 

3K 



47r oTq 



3 



For an energy, E, in the neighborhood of 10^^ ergs, a number of shells, N^, 
of 20 and a radius, tq, of 10^ cm, the temperature is below 10 MeV. 

When a shell starts, the energy Eg is very upper to the baryon mass 
energy. We can define the ratio, rj, between these two energies which is 

- _ E 

~ {M,/N,)c^ ~ ^ 

where Mf, is the total baryonic mass ejected. 

In the observer frame, the shell thickness, Ar, is su pposed to remai n 



in tne oDserver irame, tne sneii tnicKness, zir, is su pposea to remai n 
constant and equal to tq until the broadening radius (jOoodman ( 1986t ). 



Meszaros et al. ( 1993h ). Beyond this radius, defining the Lorentz factor, F, 



of the baryonic matter, the thickness becomes 

Ar ~ (r/c) Aw ~ r/2F^ (1) 

In the same frame, we also define a radius, r^, where the kinetic energy 
of baryonic matter reaches its saturation value. At this moment, the Lorentz 
factor F = T^ax is close to r]. 

In the co-moving frame, according to the Lorentz transformation, the 
shell thickness is given by 

Ai? = F Ar . 

If we consider an adiabatic expansion of the shell, conservation of entropy 
in the co-moving frame for a radius lower to r^ is given by 

r^ F ro = constant . (2) 

At the same time, conservation of energy is such that 

r^ F^ ro T'^ = constant . (3) 

From Q and Q, we deduce two laws of evolution for r < r^ which are 
F(r) oc r and T{r) oc r^^. Then, the saturation radius r^ can be defined by 
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For 7] of the order of 300, we obtain r,, ~ 3 x 10^ cm. 

Beyond the radius r^, the equation (j2I) is always vahd but the Lorentz 
factor of the shell remains constant. In this case, the temperature is such 
that T(r) oc r~^/^. Considering Ar > tq, the radius r^, according to (Q), 
satisfies Vf, > V^vq. At last, T is supposed to reach its saturation value, r], 
around well before ri,. In fact, rj is the average value of the saturation 
bulk Lorentz factor, and we have to bear in mind that deviations from this 
average is expected to generate internal shocks. The radius is so defined 



For 7] of the order of 300, we obtain ~ 9 x 10^^ cm. 

To conclude, all the different parameters of a shell we need, can be sum- 
marized by the next following expressions : 

In the observer frame, the thickness of a shell is 



by 




max 



r/rf 



for r < Tb 
for r > Tb 



The Lorentz factor of a shell is such that 




for r < Ts 
for r > rg 



In the co-moving frame, the thickness of a shell becomes equal to 



AR 



r To — r 
^max ro^r]ro 
r/Tmax ^r/r] 



for r < Tg 

for Ts < r < Tb 

for r > rb 



(4) 



The temperature of the pair electron-positron plasma is such that 



T(r) 




for r < r 



for r > r 



s 



s 
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and the co- volume of a shell, Vc{r), is 



fl r for r < r 



s 



Vc{r) ~ nr]ror'^ for Vs < r < , (5) 

(fi/r/) for r > ri, 

where Q (~ 47r/500) is the opening angle of the emission. 

The collision date (if any. . . ), tc, between two shells of Lorentz factor 71 
and 72, which started at two different times separated by Ato is such that 

2 2 

tc — 2^^^^Ato, which leads to collisions at various distances Vc such that 



rc/fb ~ Ato /to (see iDaigne fc MochkovitchI (|l998l )). For a long GRB having 



a duration ty^ ~ 10 s, the maximum collision radius reaches the deceleration 
radius r^. The light curve duration could be due to either the duration of 
the wind or to the time spread of the shell at the most remote collision, 
namely, Atd = td/iif- It turns out that this latter time is comparable to 
the duration of the wind of a long GRB, namely, 1 — 10 s. The observed 
millisecond variations should come from internal collisions located at a few 
Tfe where Ar^ ~ Tq. 

Anyway after a few r^, all the shells mix up and form a single jet. That is 
the reason why volume and thickness of a shell can really be defined only 
before the broadening radius. After this radius, in the equations (jH) and (0), 
they are undervalued. 



2.2 Conversion of the GRB energy into cosmic rays 

During the evolution, the energy of the fireball is shared between several 
forms : the thermal energy Eth-, the magnetic energy Em-, the bulk kinetic 
energy, the cosmic ray energy E^, and the energy in the form of hydromag- 
netic perturbations E^. The energy share E^ is the reservoir for particle 
acceleration. 



Until the saturation radius r^, the ratio of the pair thermal energy over 
the fireball energy, (aT^V^)/i?s, starts with a value close to unity and then 
decays as 1/r. In this stage, the proton population is a tiny contribution both 
in number of particles and in energy. Assuming that the protons dominate 
the baryonic load, we get the total number of protons 

— ~ 0.67 X lO^M — 1 



^ nip V3OO/ VlO^^er^ 
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Thus, the density of protons by shell, np{r), in the co-moving frame and 
before the broadening radius, is such that 



np{r) > 10-^ I '^^ \ I '^^ \ 



1ms J \10s J Vc{r) ' 



and the maximum value of np{r) is equal to Np/Vc{r). 
For to — tq/c — 1 ms, we can write 



3 

3 



2.7xl0-(^)-^(3^)-^(^)" (to*^)(-)" cm 

1.0 X 10- (^)-^ (3^)-^ l^y' (^) ( -^r%m-3 



np{r) > 



respectively for r < and < r < r^. 

Considering 77 = 300, even if the GRB has a long duration about 10 s, the 
density of protons, at r = r^, is larger than 10^ cm~^. We note that, up to 
the saturation radius, the density of a shell is comparable to the density of 
a sohd. 



The corresponding ratio of the thermal energy of protons, 3T / [27] mpC^) , 
starts with a very low value (~ 10~^) and decays hke 1/r as well. The 
cosmic ray component is supposed to develope out of the proton thermal 
component and would grow up to several percents of E to account for a 
significant contribution of the UHE cosmic ray generation in the Universe. 
Assume that a fraction of the proton population is injected into the cosmic 
ray component (i.e. — C*Np) and that they reach a mean energy e(r) = 
^{r)mp(? in the co-moving frame. Then, the contribution to the energy of 
the fireball, in the observer frame, is given by 

1, 

When the radius r reaches the saturation radius r^, the energy of the cos- 
mic ray component is equal to C*7(r) E/{1 + 4^*7(r)/3). For a cosmic ray 
spectrum in e~^, 7 = \og{'^rnax): which clearly shows that the goal of con- 
verting about 10 percents of the fireball energy into cosmic rays energy is 
achieved when a sizeable fraction of the protons are injected in the cosmic 
ray population. 

Regarding the magnetic energy of the fireball, two points of view can be 
considered; either the magnetic field behaves like in a jet and the poloidal 
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component decreases like thus its pressure decreases like the relativistic 
thermal pressure in whereas the toroidal component decreases more 

slowly; or the magnetic field of the shell disconnects from the central source; 
which is very likely. In this latter case, the magnetic energy of the shell 
is conserved as long as the dissipation is negligible (e.g. B oc Vc ) and 
decays when the Fermi processes and/or reconnections become efficient in 
accelerating particles. Initially, the intensity of magnetic field is usually 
considered as of the order of the equipartition value, which means that the 
magnetic energy is a sizeable fraction of the fireball energy (i?m(ro) = imE). 



2.3 The importance of pp-collisions during the primeval 
stage 

We begin with defining the radius where a shell becomes optically thin 
with respect to Compton scattering. It can easily be checked that a typical 
shell width Ai? becomes smaller than the fiow transverse radius after a short 
while, when r > rj^Jn / AVItq which is comparable to r^. It will turn out that 
the photosphere is located at a much larger distance for large enough 77 (see 
©) and therefore the opacity of a shell is determined by its width. The 
Compton opacity is = ay rie AR where ctt is the Thompson cross section. 
We assume — Up c::^ rib and thus the co-moving baryon density is related 
to the wind mass fiux by 

= FrT^ 

M Z rUp c 

Moreover for r > r^, the GRB energy is converted into relativistic kinetic 
energy such that the kinetic luminosity (kinetic energy fiux) of the wind 
= E/t^ = r] c^. For r > r^, F ~ 17 and the Compton opacity is thus 

(Tt Lyj AR 

^ 

Q nip rj"^ 

We can define a critical value for t] such that the photospheric radius is 
located at r^,, where shock acceleration starts. This critical value rj^ is given 
by 

^* \nmpC^roJ ^yiO^^ergJ \lsj V2 x IQ-V 
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The photospheric radius is such that 




for 7] < r)^ 
for Tj >rji, 



(6) 



These simple formulae illustrate the requirement of a large value of the 
asymptotic bulk Lorentz factor r] in order to observe an optically thin X-ray 
spectrum. For vj larger than the critical value, the internal shocks generate 
only a non-thermal spectrum just after r;,, whereas, for r] < r)^, the internal 
shocks start accelerating particles in an optically thick plasma. Assuming 
that the GRB energy varies within two orders of magnitude, the wind time 
by two orders of magnitude and the solid angle by one order of magnitude, 
the possible values of rj^ extend within an interval of one order of magnitude. 
However it is worth mentioning that this usual approach of the opacity issue 
is an underestimate of the opacity effects, because the coalescence of shells 
and ultimately their dissolution into a single jet increase the opacity. 

The possibility for the internal shock scenario to partially operate in an 
optically thick regime is important when one consider pp-coUisions. Indeed, 
any relativistic protons can collide with other protons (or neutrons) produc- 
ing pions and thus neutrinos. The cross section, CTpp, for this production is 
constant and equal to 2.7 x 10^^^ cm^ when the kinetic energy of protons is 
larger than 1 GeV. The sheets become thin to pp-coUisions always before the 
photosphere since Tpp = Hp app AR = 1 at Vpp such that Vpp — r-^^fop^foT < r^, 
for Tpp < Tft and r^p = r^app/ax < for rpp > r^. In this latter case, 
Tpp ~ 0.04 r^; which occurs when 1] < ?7^/2. In this paper, we intend to em- 
phasize the importance of this energy limitation to proton acceleration and 
estimate the resulting neutrino emission. 

3 Consequences of a fast acceleration regime 

In this section, we analyse the consequences of the usual assumption of a 
Bohm scaling for the Fermi process for both the energy distribution cut off 
and the depletion of the energy reservoir for acceleration. 
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3.1 Bohm scaling and energy losses 

We first consider a Fermi acceleration process that is assumed to follow the 
so called Bohm's scaling namely characterized by its time scale proportional 
to the Larmor time, i.e. : 

where we take Kq = 10, as often assumed in astroparticle physics. 

In the early phase, the cut off of the proton energy distribution can be 
caused by cither the synchrotron loss or the pp- collisions; which is obtained 
by equating the acceleration time and the loss time. The synchrotron time 
tsyn is given by 



(^) w^, 



1 _ AZ^aT^ (m, 
where Wm is the density of magnetic energy. 
This expression leads to the cut off energy : 

1 f B\ 

e,j,„~2.4xl0i^^ — GeV . 

We consider now the limitation caused by pp-coUisions. The cut off energy 
Cpp is such that {upappc) Kq^l = 1) and thus 



At the broadening radius (rt), where the magnetic field could be as high as 
10^ G, these two limitations are comparable : 

and 



107 g; V300/ \n. 
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During the expansion the ratio between these two hmitations evolves such 
that 

^ 1.2 X 10^ m ( 

esyn VIO/ \300J \riJ yiQ-'G J \rj 

which suggests that for any decrease of the magnetic field in r~" with a 
between 1 and 2, the most severe limitation would be due to synchrotron 
losses. However we will show in the next subsection that such efficient Fermi 
acceleration would dissipate the magnetic energy very rapidly. 



3.2 The depletion of the energy reservoir for acceler- 
ation 

Such an efficient acceleration produces a strong depletion of the energy 
reservoir for particle acceleration. The previous results would make sense 
only if the depletion of the energy reservoir for particle acceleration is slow. 
This is exactly the purpose of this subsection. 

Let be the energy in the hydromagnetic perturbations (magnetic and 
kinetic energy) involved in the acceleration of particles. We assume that 
this is a sizeable fraction of the total magnetic energy and that even the 
large scale magnetic field in the shell is mostly disconnected from the central 
source after a short while and thus is rather tangle. Therefore we assume that 
the total magnetic energy dissipates at the same rate than E^. A minimum 
dissipation is assumed by considering its energy loss by particle acceleration 
only. According to (13. IB . the acceleration power is = -^^m„c^, where 

t^Q is the Larmor time for a proton of 1 GeV in the co-moving frame. The 
minimum depletion is then governed by E'^ = —Q^, which reads : 

E B 

^™ = -^^, (7) 

where we define a characteristic time t'^ at a given radius f where the magnetic 
field has an intensity -B^,, such that 

Clearly, this time measures the rate of depletion at the considered radius f 
and must be compared to the co-moving dynamical time t'{f) = r/Tc. Let us 
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make an estimate at r;,, the late stage actually; for B = 10^ G and e = 1 GeV, 
t'^ ~ 3x 10~® s, whereas t'{rh) ~ 10~^ s. This perturbation burning out within 
a few nanoseconds is not realistic at all. Before this stage, the depletion time 
is much shorter than the fireball dynamical time by a factor proportional to 
(Br)~^ that goes at best like r^/^ for B oc r'^^"^ corresponding to magnetic 
energy conservation and the perturbations are burnt out almost immediately, 
their energy being radiated by photons and neutrinos. The observation of 
these perturbations in the light curve of gamma emissions suggests a much 
slower depletion; which we will analyse in the next section. 

4 The interest of a progressive acceleration 
regime 

The Bohm scaling is often used in astroparticle physics, which sometimes 
can pro vide some rough e stimate of the high energy cut off. But the true 



scaling (|Casse et al.l 120011 ) depends on the turbulence spectrum. We will 
assume that the perturbations are distributed according to the Kolmogorov 
law. This will significantly change the performance of the Fermi process and 
the conclusions about cosmic ray generation. 

4.1 Kolmogorov scaling and energy losses 

In the co-moving frame, the acceleration time depends on the speed c 
of magnetic perturbations that scatter particles on both sides of the shocks. 
Considering the same average velocity, we have the relation 

ts 

f = — 

f'acc ^2 ' 

where the scattering time, tg, can be expressed like 

with rjb = ^^§2^ and p = < 1 where Ic is the mean coherence length of 
the magnetic perturbations whose spectrum is supposed to be a power law 
of index j3. According to the Kolmogorov theory, one can take /? = 5/3 and, 
with kq = l/{P'^r]b), we get 

K, = kq p~^/^ . 
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Considering Ic < Ai?, 



and we have 



AR ~ ZeBARc 

2/3 / A D \ 2/3 



The synchrotron hmitation becomes 



For r = Tfe (~ 10^2 cm), 5 = 10^ G and = 47r/500, e,j^„ ^ 3.8 x 10^ GeV. 



As for the pp-colhsion hmitation, 



^^oY^ ( T] \^ ( r]V^ f B \ f AR \''^fr^^ 



For r > Tfe, AR = r /rj, defining r?,^ = ?7^ro, we obtain 

This energy increases up to a maximum value reached at r^p (when it's larger 
than Tf,), namely. 

As long as r < r^p, the energy limitation is clearly due to pp-collisions rather 
than synchrotron losses. These results show, first, that it is not possible 
to accelerate protons beyond 100 GeV before the fireball becomes thin for 
protons at r^p, for 77 > 77^; second, that, for rj < ri^jl^ the proton energy could 
increase above 100 GeV. For instance, if 77 = ?7*/3, ~ 240 r^, rpp ~ 10 r^, 
B weakens by a factor {rpp/rf,)'^ ^ 3.2 x 10~^ with a = 3/2 and thus 
compensates the increase due to rj. Therefore the energy remains of order 100 
GeV. We have to examine how this regime changes the rate of perturbation 
energy depletion. 

In the opposite situation where the protons are not accelerated beyond GeV 
in the opaque stage, for r] > rj^ , and /or < 10^ G at r^, nucleosynthesis is 
possible as shown by Lemoine ( Lemoinel 2002| ) . 
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4.2 The depletion of the energy reservoir for acceler- 
ation 



Taking into account the variation of k in the acceleration power, we define 
a new characteristic time t'^ which is quite different from the last one : 



For a power law distribution, not harder than e~^, < p^/^ >~ Pq , po being 
the rigidity for 1 GeV. For bell type distribution of standard deviation e, 
< p^/^ >~< p(e)^/^ >. Compared to the Bohm regime, the Kolmogorov 
regime increases this time by a factor Pq . At Vb we get a more extended 
time of the order of 0.1 s, which is much more realistic than in the Bohm 
scaling case. If gnificantly larger, the depletion behaviour is changed 

such that t^/t' decreases instead of increasing like in Bohm regime. Indeed in 
Kolmogorov regime, this ratio scales like B~^^^r~^^^, which decreases when 
B decreases less slowly than r~^/^. 

For a more detailed estimate, we integrated the differential equation d?)) 
which, after some algebra, leads to the following result : 



In fact, the perturbations are burnt out in a fraction 6t' of t'^{rb), which can 
be obtained directly from the differential equation ^ : 



which confirms the previous statement of a reasonable depletion time. 

4.3 The pp-neutrino emission 

In this subsection, we give an evaluation of the number of emitted pp- 
neutrinos and we calculate their energy spectrum. 




St' = t'M 



E 
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The pp-collisions produce neutrinos after some reactions which are 

P + P > + TT"*" 

p + p + a(7r+ + 7r^) + hn^ 
p + n + vr'*' + a(7r^ + tt^) + Im^ 
2n + 27r+ + a(7r+ + tt") + 67r° 

The production of mesons vr" and tt"*" gives neutrinos et /x-mesons through 
the decay reactions 



TT > jj, + 



TT 



+ 



/i+ + Z/^ (11) 



and other neutrinos are produced after 



/i+ — ^e+ + z/e + z7^. (12) 

Let's write energy of different particles 

m^± ~ 140 MeV/c^ 

- 105 MeV/c^ . 
me± ~ Q.^MeV/c^ 

Thus, neutrinos coming from have a minimum energy equal to 

2 2 

m ± — mf, 
~ ^ ^ ~ 30 MeV 

2 m^± 

and decay of yU-mesons (II 2j) gives neutrinos having a minimum energy which 
varies from 25 to 50 MeV. 

Let fv{^v) be the energy distribution of neutrinos, normalised such that 
/ fui^u) dey = 111,, the number density of neutrinos; and /*(7) the Lorentz fac- 
tor distribution of relativistic protons, normalised such that / /*(7) d'j = n-^, 
the number density of relativistic protons. The kinetic equation for neutrinos 
can be written in a simplified way as follows : 

d f 

—fy + cn.Vfy = Vpp j f^i'j) 5{e^ - 7ceo) dj , (13) 
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where Sq is the average energy of the neutrinos generated by the pp-colhsion 
in the center of mass frame (above estimated), the average number of 
produced neutrinos at each coUision and 7c the Lorentz factor of the coUision 
frame. The energy of the emitted neutrinos remain close to the minimum 
values because the pions do not take all the energy available above the thresh- 
old, but just a little excess above it, the remaining energy being kept by the 
proton or the neutron. Therefore the number is just a few. The Lorentz 
factor 7c ~ ^/7, where 7 is the Lorentz factor of the relativistic proton. So 
we can easily deduce that, for oc 7"'', the energy flux of neutrinos is in 
e"^^"*"^, between eo and eQ^fe^^^fm^ ^ in the co- moving frame. For s = 2, a 
neutrino spectrum in is expected between 5 GeV and 150 GeV, typically, 
for the observer frame seing the GRB flow coming front. 

Assuming an isotropic neutrino emission in the co-moving frame, the 
power emitted is 



where the integral over the width extends down to the lower bound at 
Tpp, corresponding to Tpp = 1. Thus E'^, = t'^Vt^y eQcni^r'^^. The quan- 
tity t'^Q cn^^r"^,^ is the number of relativistic protons ^i,Np which have flown 
through the "proto" -sphere during the GRB event. Finally, we get a simple 
formula giving the amount of neutrino energy in the co-moving frame : 

El^^.C.^-. (15) 

For an observer seing the GRB shells coming front, the neutrino energy 
emission is multiplied by the bulk Lorentz factor, thus = ,^^^i^(eo/mpC^)£' 
which reasonably leads to E^, ~ lO^'^ — lO^^E. The number of emitted neu- 
trinos by a GRB is simply ^^^-k- A neutrino telescope of collecting surface 
A can receive N^A/AttD'^ from a GRB exploding at a distance D. For ex- 
ample, a GRB occuring at 100 Mpc can provide with 10^ neutrinos crossing 
a A:m^-dctcctor. However, the number of events recorded by the detector is 
obtained by multiplying this number by the detection probability which is 
quite low. 

In the same process, vr'^-decay occurs and generates gamma photons. An 
excess of these photons could be seen in the gamma spectrum; they superim- 
pose a spectrum in e;^;^*'^^ on the synchrotron spectrum due to the electrons. 




(14) 
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Because we focused on proton acceleration and neutrino emission, we disre- 
garded the interesting issue of the neutron component which can decouple 
under some conditions ( Derishev et al. 19991 ). 



The stage of UHE Cosmic Rays accelera- 
tion 

Because of their magnetic field intensity and the size of their wi nd, the 



GRBs are consid e red a s possible accelerators of UHE-cosmic rays (|Vietri 



f 19951 ). Waxman ( 19951 )) : for r > r^, particles could reach the maximal 



energy beyond which they are no longer confined, emax = Z eB r, namely 

emax = 7.7 X 10^1 (t^) (^T^W-z:) eV^> 
for a proton [Z = 1). 



10^ G V 1012 cm 



5.1 Proton energy in the transparent stage 

We examine the loss limitation suffered by the protons after crossing the 
"proto" -sphere. After the equation (jH)), we have already given the value of 
the maximum energy due to the synchrotron loss at r^,; thus, it scales like : 



5/4 / ^ ^ (5q-2)/4 



It will turn out that this value is above the threshold of the p7-process, but 
is far below the range expected for getting UHE-Cosmic Rays. For a = 3/2, 
this range would be reached at a few 100 ri,. 

The expansion loss can be significant. Indeed, the cut off energy is ob- 
tained by setting tacc = t' = t/T in the co-moving frame, and the limitation, 
eexp, is such that 



e ~ 



Whereas the synchrotron loss diminishes with distance, the expansion limi- 
tation becomes more and more severe (see figure Q). Moreover, it can easily 
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be checked that the escape time due to transverse diffusion is always much 
longer than the expansion time. 




1 10 100 r/x^ 



B=10^x(r/rJ G 

b 

Figure 1: Diagram of energy limitation in the co-moving frame for a = 3/2 : 
the dashed lines show the result of the Bohm scaling and the solid lines are 
obtained with the Kolmogorov scaling. The horizontal dashed line represents 
the lowest UHECR range limit. 



5.2 The p7-neutrino emission 

A signature of VHE-protons acceleration is expected with the neutrino 
production resulting from the photo-production of pions, which is efficient to 
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produce neutrinos through the A — resonance (the so-called GZK effect) 



(2/3) n 

P + TT + 7 + 7 

p + 7 (18) 
n + 77^ ^ ... ^ p + e^ + e + + J^e + + J^iJ. 

The threshold of the process is such that e-y > \mnniT, — jml, which leads 
to a threshold for the proton Lorentz factor ^th ~ 10^, the energy of the 
photons being in the keV range in the co-moving frame. As previously seen, 
7max is at least ten times larger than 7j/j. The GRBs are likely opaque to 
this process until the shells reach the radius r^^ such that 

y P X 1/2 

where e-y is the average energy of the target photons in the observer frame 
(typically 1 MeV); which puts this p7-"proto" -sphere at a rather large radius 
of 1.5 X 10^^ cm for = 0.1 E. 

In the collision frame, the energy of a pion, generated sufficiently above 
the threshold, is ~ ~ el^, ~ 767(1 — cos 6'), where 7 is the Lorentz factor 
of the proton. The emitted neutrinos have an energy which is a fraction 
ao (~ 5%) of the pion energy. Similarly to the statistical treatment of the 
pp-collisions, we write a simplified (with delta approximation instead of a 
function smoothed by angle averaging) kinetic equation for the isotropic dis- 
tribution of the neutrinos generated by p7-collisions in the co-moving frame, 
which is : 

+ cn.Vf,y = ^i^cajyy j /*(7) ^7 / /^(e^) 5{e„ - a^^'^e^) de^ . (20) 

Integrating both sides of the equation over the co-volume and the proper 
time, we easily get the number of emitted neutrinos : 

N,^i,N> ^i^nj- ^"j , (21) 

where is the number of cosmic rays above the p7-threshold energy, the 
last result being obtained for a 7"^ proton distribution. The integration, like 
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in the case of pp-collisions, starts at the p7-"proto" -sphere. The total energy 
radiated by neutrinos in the co-moving flow is also derived easily from the 
kinetic equation : 

El-aoWog(^) 6-;. (22) 
\ 1th J 

The logarithm factor is obtained for a distribution of protons and can 
easily be modified. This result can also be rewritten in the observer frame 
as follows : 

Ey ~ ao ith £7 log ( ) Ny . (23) 
V 1th J 

The energy radiated in the form of p7-neutrinos is E,j ^ 10"'' - E. Nu merical 



computation of such spectra has been done (see iMiicke et al. (|l99i l. The 
kinetic equation also provides with the neutrino spectrum. Assuming is a 
power law distribution and that the energy distribution of target photons is 
oc 6"°", we get a power law energy spectrum for the neutrinos, namely 



2 dNy{> ey) 



oc tu ^ 



in an energy range depending on the energy range of the protons since e,^ ~ 
aoi'^e^. These spectra are given in the co-moving frame. For the observer, 
they are Doppler beamed with the bulk Lorentz factor of the relativistic 
wind. However, the number of events is so low that it is still "virtual" to 
talk about a spectrum... 



6 Discussion 

In order to account for the non-thermal and highly variable gamma emis- 
sion of GRBs, the fireball model and the internal shock model have been de- 
signed with baryonic load parameter rj that has been supposed large enough 
to get a relativistic wind achieving a large Lorentz factor F ~ 77. However, it 
turns out that opacity effects could easily be significant at the beginning of 
the emission. We stress that point through a discussion involving a critical 
value rji, of the baryon load parameter. Indeed one emphasis of this paper 
is to analyse the opacity of the GRB to relativistic protons with respect to 
pion production by pp-collisions. Actually, the opacity condition relative to 
pp-collision is not far from the Compton opacity condition. Therefore it is 
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reasonable to think that a significant fraction of the GRBs experience proton 
Fermi acceleration with efficient pp-coUisions revealed by neutrino radiation. 

During a stage of pp-opacity, we have shown that Fermi acceleration leads 
to conclusions about the pp-process that are very sensitive to the choice of 
the efficiency law of the acceleration. We have shown that the Bohm scaling 
assumption leads to an efficient proton acceleration that would be limited by 
synchroton loss, whereas the correct law governed by the turbulence spec- 
trum (Kolmogorov law was used) leads to very different estimates and the 
pp-collisions process turns out to be the main limitation of the proton en- 
ergy. Moreover, the excessive efficiency of the Bohm scaling would make the 
acceleration to deplete its energy reservoir in a time too short to maintain 
the gamma emission. The most important surprise raised in estimating the 
energy limitation in the radiation free stage. Indeed, the expansion limitation 
turns out to be drastic with the Kolmogorov law and maintains the proton 
energy below a few 10^ GeV. 

The analysis of the paper indicates that a double neutrino emission can 
be expected with many GRBs, namely a stage of pp-neutrino emission fol- 
lowed by a stage of p7-neutrinos. The number of emitted pp-neutrinos gives 
the amount of relativistic protons; and the number of p7-neutrinos gives the 
number of protons above the threshold, which is about 10~^ less. We pro- 
posed an analytical shape of the neutrinos spectra, as well as the photon 
spectrum generated by the 7r°-decay. 

The magnetic field intensity is an important parameter that controls the 
proton acceleration and the synchrotron losses. The neutrino emissions sig- 
nificatively depend on its value at the crucial distance Vb- We took a high 
but still reasonable value of lO'^G at this distance; if we take less, the syn- 
chrotron limitation is less important, but the expansion limitation, which 
controls the highest energy of the protons in the GRB, becomes more se- 
vere, and also the pp-neutrino emission becomes less energetic. If we unduly 
take more, synchrotron losses dominate over pp-collisions and the accelera- 
tion is more efficient against the expansion losses, however not sufficient to 
get UHE-cosmic rays. . . This paper does not exclude the possibility of UHE- 
cosmic ray generation in GRBs. It simply states that its achievement with 
Bohm scaling is not reliable and leads to observational inconsistencies and 
that its achievement with Kolmogorov scal ii ig is impossible. . . We thin k that 
there is another possibility (|Pelletiei (|l999l ). IPelletier fc Kersali (|2nnnl )) that 
deserves a detailed investigation that we will present in a forthcoming paper. 
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